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Abstract The nucleus of the Seyfert galaxy NGC 1068 is believed to host a supermassive black hole. Evidence for 
the presence of a massive central object is provided by water maser emission, which displays a linear pattern in 
the sky, suggestive of a rotating disk. The rotating disk hypothesis is further strengthened by the declining shape 
of the derived rotation curve. Similar maser emission from NGC 4258 has led to a reliable estimate of the mass of 
the central black hole, because in this case the rotation curve is Keplerian. In the case of NGC 1068 the rotation 
curve traced by the water maser is non-Keplerian. In this paper we provide an interpretation of the non-Keplerian 
rotation in NGC 1068 by means of a self-gravitating accretion disk model. We obtain a good fit to the available 
data and derive a black hole mass M, ~ (8.0 ± 0.3) 10® Mq. The resulting disk mass is comparable to the black 
hole mass. As an interesting by-product of our fitting procedure, we are able to estimate the viscosity parameter, 
which turns out to be a ~ 10“^, in line with some theoretical expectations. 
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1. Introduction 

There is now convincing evidence that most AGNs host 
a supermassive black hole, with masses ranging from 
IO^Mq to 10® Mq. The gravitational energy extracted 
from an accretion disk around such black holes is generally 
considered to be the main source of the AGN luminosity. 
Determining the mass of the central black hole M, is thus 
one important goal of studies of AGlNs. This has received 
greater attention recently, in a more general context, be- 
cause correlations have been found between M, and the 
global properties of the host galaxy, such as its mass or 


luminosity (Magorrian et al., 

1998 

) or its central veloc 

ity dispersion (Ferrarese & Merritt, 

200C; Gebhardt et al. 


2000a). 


The central point mass M, can be measured by differ¬ 
ent methods. The analysis of JIST optical spectra in terms 
of stellar dynamical models (see, for example, |Gebhard~ 


sight. An independent powerful tool to obtain reliable esti¬ 
mates of M, is provided by the study of the Doppler shift 
of water maser emission lines (when available). This latter 
method provides the most reliable determinations of M,, 
because it probes the gravitational held at very small dis¬ 
tances (less than 1 pc) from the center. The masing spots 
are often observed to trace a linear structure and show a 


declining rotation curve. In the case of NGC 4258 (Miyoshi 


et ah, 1995), the rotation curve is remarkably Keplerian 
and leads to a robust determination of M, k, 3.9 IO^Mq. 

NGC 1068 is one of the best studied Seyfert galaxies. 
It is considered to be the prototypical Type II Seyfert, in 
which the central engine is believed to be hidden from our 
line of sight by a dusty structure. This obscuring structure, 
now resolved with VLBA radio continuum observations 


et ah, 2000b) and the study of gas kinematics from HST 
spectra (van der Marel & van den Bosch, 1998|) probe the 


nuclear gravitational held of nearby galaxies at distances 
typically of the order of « 100 pc from the center. M, 
can be estimated also from reverberation mapping of the 
broad line region, which is able to probe the gravitational 
held at smaller scales. However, reverberation mapping is 


( Gallimore et ah , 1997 ), appears to have a disk-like rather 
than a toroidal shape. The bolometric luminosity of the 
nucleus is Lboi ~ 8 10^"^ erg/sec ( Pier et ah , 1994 ). If 
we adopt for the accretion efficiency a value of ry = 0.06, 
appropriate for a non-rotating black hole, the resulting 
mass accretion rate is M « 0.235 Mq/ yr. Throughout 
this paper we will assume that the distance to NGC 1068 
is 14.4 Mpc (so that 1 mas = 0.069 pc). 

Water maser emission is observed from the nucleus of 
NGC 1068 extending out to a radius r « 1 pc (Greenhill 


not applicable to Type II Seyfert galaxies, like NGC 1068, et ah, 1996; Greenhill & Gwinn, 1997). This emission is 


in which the broad line region is hidden from our line of believed to come from a rotating, almost edge-on accre- 
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tion disk, possibly associated with the disk seen in VLBA 
continuum (Gallimore et ah, 1997). The striking feature 
of this maser emission is that, in contrast with the case 
of NGC 4258, the rotation curve of the masing spots is 
clearly non-Keplerian. preenhill et al. ( 1996| ) report a best 
fit to the data with a circular velocity V oc 7--o.3i^ Possible 
causes for the sub-Keplerian rotation may be: a) the maser 
may not be associated with a pure rotating structure, be¬ 
cause an outflow is present; b) radiation pressure may re¬ 
duce the rotational velocity (Pier & Krolik, 1992a); c) the 
source of the gravitational field may be extended, includ- 


fies the vertical hydrostatic equilibrium (Paczyhski, 1978; 


2.1. Self-regulation of Jeans instability 

The onset of gravitational instabilities in a fluid disk is 
controlled by the well-known axisymmetric stability pa¬ 
rameter Q: 


Q = 


CsK 

irGa ’ 


( 1 ) 


ing a nuclear stellar cluster (as proposed by Kumar 1999) 
or the accretion disk itself. 

In this paper we interpret the non-Keplerian rotation 
in NGC 1068 in terms of a self-gravitating accretion disk 
model (Bertin & Lodato 199E, hereafter BL). The disk 
is assumed to be self-regulated at the threshold of Jeans 
instability. The gravitational field is computed by solving 
the relevant Poisson equation including both the central 
point-like object and the disk. The resulting rotation curve 
is characterized by an inner Keplerian curve connected 
with an outer asymptotically flat rotation curve. 

We have fitt ed the Very Long Baseline Interferometry 
(VLBI) data of Greenhill fc Gwinc ( 1997| ) with our the¬ 
oretical models. The fit is very good. From the results of 
the fit we can derive the value of the central black hole 
mass, M, ft! 8 10®Mq, roughly one half of the value that 
would be inferred from a Keplerian fit (from which one 
would obtain M, ft 1.5 IO^Mq), basically because we at¬ 
tribute part of the gravitating mass to the disk. As an 
interesting by-product of our modeling procedure, we are 
able to derive the value of the viscosity parameter a that 
regulates the accretion process, which turns out to be in 
agreement with some theoretical expectations. 

The paper is organized as follows. In Section || we de¬ 
scribe the role played by the disk self-gravity in the outer 
parts of AGNs. In Section || we report the observations of 
the nucleus of NGC 1068, focusing on its geometry and 
kinematics. In Section 0 we describe our fit to the water 
maser data. In Section |[we discuss some possible alterna¬ 
tives to the model presented in this paper. In Section ^we 
draw our conclusions. In Appendix ^ we describe in some 
detail the statistical significance of the fitting procedure 
adopted. 


2. Preliminary considerations on the influence of 
the disk self-gravity in AGN accretion disks 

The disk self-gravity may affect several aspects of the dy¬ 
namics of accretion disks: a) gravitational instabilities are 
expected to modify the energy and angular momentum 
transport in the disk, perhaps being the main tool able to 
drive accretion at large radii ( Lin fc Pringle^ 1987); b) the 
vertical gravitational field associated with the disk modi¬ 


Bardou et al, 1998); c) the radial gravitational field of the 


where Cg is the thermal speed, k is the epicyclic frequency, 
and (T is the disk surface density. For Q < 1 the disk is 
unstable. As Q is proportional to the thermal speed, a hot 
enough disk is expected to be free from the effects related 
to the disk self-gravity, while, on the other hand, if the 
disk is cold enough to begin with, it cannot survive long 
in such a condition and is likely to be eventually char¬ 
acterized by a value of Q close to unity, as a result of a 
self-regulation mechanism, studied and recognized espe¬ 


cially in the field of galactic dynamics (see Bertin fc Lin 


1996 ). Self-regulation has been sometimes taken into ac¬ 


count also in the context of accretion disks ( Lin fc Pringle 
1987; Hure, 200C|). The inner parts of the accretion disk in 


AGNs are very hot, so that the disk self-gravity is not ex¬ 
pected to play a role there. On the other hand, it is easy to 
show that the outer colder parts of the disk may be subject 
to gravitational instabilities. Consider the outer region of 
the a-disk solution by Shakura fc Sunyaev ( 1973 ), for the 
case in which gas pressure and free-free absorption domi¬ 
nate. The radial dependence of Q in this case is: 


Q ft 5.6 


( 2 ) 


disk may lead to deviations from Keplerian rotation ( 


where we have scaled the basic physical parameters to 
typical AGN values, so that M 26 is the mass accretion 
rate in units of 10^® g/sec ft 1.57MQ/yr, Mg is the black 
hole mass in units of 10® Mq, and f is the radius in units of 
the Schwarzschild radius of the black hole R, — 2GM,/(?\ 
a is the dimensionless viscosity parameter entering the 
Shakura-Sunyaev prescription for viscosity. The value of Q 
thus decreases rapidly with increasing radius and becomes 
equal to unity at fq = rq/R, ft 10® ( for Mg = M 26 = 1 
and a — 0.01), i.e. at rq ft 10“^ pc. Kumar ( 1999 ) has 
carried out a similar analysis in the case in which the 
opacity is dominated by metal grains and has basically 
confirmed the above estimate, finding that Q ft 1 at rg ft 
3 10“® pc, for the same input parameters. 

The argument of self-regulation suggests that the outer 
disk, beyond rq, be characterized by Q ft 1. 

2.2. Impact of the disk self-gravity on the vertical 
structure 

If we consider the modifications of the vertical hydro¬ 
static equilibrium by the disk self-gravity, a simple way 
to address the problem is to compare the vertical length- 
scales derived in the limiting cases of non-self-gravitating 
disk and of fully self-gravitating disk. In the non-self- 
gravitating case we have: 


hnsg - ^ ■. 


( 3 ) 
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while in the self-gravitating case: 
irGa' 

The disk self-gravity will then modify the vertical struc¬ 
ture of the disk if: 


(4) 




ttGct 




(5) 


because k is related to by a numerical factor (dependent 
on the rotation curve) close to unity. Therefore, a discus¬ 
sion of the vertical equilibrium also leads to the conclusion 
that an accretion disk is expected to be self-gravitating if 
it extends to radii larger than rg. This simplified analysis 
is confirmed by the more detailed study presented earlier 
by us (BL, Appendix), where a convenient analytical ex¬ 


ttGct 4g2(2J72/^2 _ X) 

/l +- 1 . (6) 




- 2/3 


(7) 


NGC 1068 traces the outer parts of the nuclear accretion 
disk, such disk is likely to be self-gravitating for what con¬ 
cerns its vertical structure and transport phenomena, and 
it should display significant deviations from Keplerian ro¬ 
tation. 

3. The masing spots in NGC 1068: geometry and 
kinematics 

In this Section we will discuss the geometry of the maser 
emission in NGC 1068, for which different interpretations 
have been given. 

3.1. Disk or torus? 


pression for the disk thickness, in the case in which both 
the disk and the central object are taken into account, is 
provided: 


h = 


2.3. Modification of the rotation curve 

In this paper we w ill r efer to the self-gravitating steady- 
state disk model of BL, who considered disks that are self- 


regulated with respect to Jeans instability, so that Q = 
Q re 1. The gravitational field of the disk is computed by 
solving self-consistently the relevant Poisson equation (see 
Eq. (4) of 13). 

The basic dynamical feature of such self-regulated disk 
model is that at large radii the rotation curve V(r) departs 
from the Keplerian profile, approaching an asymptotically 
flat rotation curve. The typical lengthscale that marks the 
transition from Keplerian to self-gravitating regime is: 


Greenhill et al. ( 1996D observed the “red-shifted” water 
maser emission from NGC 1068, finding that it traces 
a linear pattern in the plane of the sky, inclined by ap¬ 
proximately 45° with respect to the direction of the ra¬ 
dio jet. The misalignment between the radio axis and 
the maser emission led to an early interpretation of the 
maser spots as arising from the upper limb of a nearly 
edge-on torus, characterized by a rather large aspect ratio 
(see the schematic representation in Fig. |]). According to 
this interpretation, we would expect that the correspond¬ 
ing “blue-shifted” emission should come from the dashed 
region shown in Fig. |I|. On the o ther hand, subsequent 
observations ( [Greenhill fc Gwinn , 1997 ) showed that the 
emission traces a linear pattern from the “red-shifted” to 
the “blue-shifted” emission, hence arguing in favor of a 
thin disk interpretation. The misalignment between the 
disk axis and the ra dio jet is not uncommon in AGNs (see 
$chmitt et al. , 20021) and may be due to a variety of phys¬ 
ical mechanisms (for example, a warp in the outer disk; 
Pringle] [l9^. 


Deviations from Keplerian rotation occur already at radii 
of the order of a fraction of rg; in fact, we can have 
d InK/d Inr = —0.4 at r = O.lrg. For r Vg, the mass 
of the disk grows linearly, and the surface density behaves 
as (T oc l/r. If we estimate the value of Cg, assuming the 
typical AGN values for the relevant parameters used ear¬ 
lier (Mg = M 26 = 1 and a = 0.01), we find that Vg is of 
the order of a few pc. Interestingly, we find that deviations 
from Keplerian rotation are expected to occur exactly at a 
distance from the center probed, in the case of NGC 1068, 
by the water maser emission. In contrast, in the case of 
NGC 4258, a similar analysis would lead to an estimated 
Tg 3> 1 pc, consistent with the Keplerian rotation curve 


observed by Miyoshi et al. (1995). 

Based on the arguments presented in this Section, we 
therefore conclude that, if the water maser emission from 


3.2. Geometry of the disk emission 

The rotation curve of the maser spots can be divided in 
two different regions: 1) at small impact parameter the 
masers show an apparently rising rotation curve; 2) start¬ 
ing from a radial distance w 0.6 pc from the center the 
rotation curve declines, following a sub-Keplerian prohle. 

The natural interpretation of the declining part of the 
rotation curve is that it arises from material that moves 
parallel to our line of sight (i.e. that lies on a disk diameter 
perpendicular to the line of sight). The best argument 
in favor of this interpretation is that maser amplification 
is largest for material that lies close to the line of the 
nodes. On the other hand, the rising part of the observed 
“rotation curve” is thought to originate from one quarter 
of the disk at the inner maser disk radius, so that the 
rising curve is an effect of velocity projection along the 
line of sight (see also Miyoshi et al., 1995). According to 
this interpretation, the inner radius of the maser disk is 
located at w 0.6 pc and the outer disk radius is at w 1 pc. 


Baan & Hashlnck (1996) claimed to have observed a 
drift in the velocity of the water masers, indicating that 
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Figure 1. Schematic representation of the location of the maser emission according to the torus interpretation. The 
three elliptical patches show the location of the observed “red-shifted” masers, of the pre dicted “blue-shifted” e mission, 
based on the torus interpretation, and of the observed “blue-shifted” emission (see also Greenhill et al. 1996). 


the maser spots are subject to large accelerations, incom¬ 
patible with the disk interpretation, according to which 
the maser centripetal acceleration should be perpendicu¬ 
lar to the line of sight. Those large accelerations have not 


been confirmed by subsequent work by Gallimore et al. 
( |2001| ), who, monitoring the velocity drift, hnd that the 
maser spots between 0.6 pc and 1 pc should lie within 
9 < 2° from the line of nodes. 


In the following we will therefore assume that the water 
maser emission comes from an edge-on thin disk extending 
from 0.6 pc to 1 pc and that the declining part of the 
rotation curve comes from material that lies within 2° 
from the line of nodes. 


4. Non-Keplerian rotation as an effect of the disk 
self-gravity 


lengthscale (defined in Eq. (|^)), the velocity scale Vq, 
defined as: 




/ • \ 2/3 

8 I GM\ 

^ ) 


( 8 ) 


and a dimensionless parameter (called ^), which is propor¬ 
tional to the net angular momentum flux in the disk J. 
We will fix this parameter by requiring the no-torque con¬ 
dition at the radius corresponding to the innermost stable 
orbit around the black hole. However, the specific choice of 
^ is not critical for our conclusions, because we are explor¬ 
ing the disk properties at large radii, where the effects of 
the inner boundary condition are negligible. Therefore, we 
are left with two free dimensional scal es, that are obtained 
by fi tting the available data (from [Greenhill fc Gwiim 
19^ . 

Note that, from the definitions of rg and Vqi we have 
that: 


4.1. Fit procedure and results 


T/2 _ 

Vq — 


GM, 


(9) 


Here we will fit the kinematical data of NGG 106 8 by using 
the completely self-regulated disk model of BL (in which, 
for simplicity, the self-regulation prescription is taken to 
hold at all radii). Actually, we expect the inner disk to 
be hotter and characterized by a higher value of the sta¬ 
bility parameter Q, so t hat a partially self-regulated disk 
model (also de scri bed in M) would be more appropriate. 
In Subsection 4.3 we will justify the consistency of the 
simpler model adopted below. 

The rotation curve of the completely self-regulated 


Therefore, from the value of Vs and Vq we can estimate M, 
and M/a. In addition, once these quantities are specified, 
from the models of BL we can obtain the disk mass. 

In our fit, we have restricted our attention to the data 
of the declining part of the rotation curve of the “red- 
shifted” maser (i.e. the data points with r > 0.6 pc), be¬ 
cause these data are directly related to the gravitational 
potential. In contrast, the velocity data points at r < 0.6, 
according to the interpretation of the masing disk geome¬ 


try presented in Subsection 3.2, only reflect the rotation at 


model by BL is determined when one specifies the radial the inner radius of the masing disk and do not carry any 
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Radius (pc) 


Figu re 2. Fit to the rot a tion c urve from the water maser emission by a self-gravitating accretion disk model. Data 
from Greenhill fc Gwinn ( 1997 ). The small panel shows a blow up of the declining part of the rotation curve together 
with the best fit obtained by assuming Keplerian rotation. The error bars reported here include uncertainties beyond 
the actual VLBI spectral resolution (see text). 


additional information about the mass distribution. We 
have assumed that the sys temic velocity is 1126 km/sec 
( [Greenhill fc Gwinn , 1997 ). T he uncertainty in th e posi- 
tion of the maser spots in the [Greenhill fc Gwiimj ( 1997 ) 
data is of the order of ~ 50/ias. The spectral resolution 
of the VLBI data is < 1 km/sec. However, the velocity 
uncertainty AV is expected to be higher due to the un¬ 
certainties in the estimate of the systemic velocity and of 
the magnitude of turbulent motion in the disk. We will 
assume that these uncertainties sum up to « 10 km/sec 
(see also discussion in Appendix . All fits are obtained 
with a y-square minimization. The result of the fit, shown 
in Fig. is satisfactory. The resulting reduced y-square is 
= 0.55, with 48 degrees of freedom. The fit parameters 
are Vq = (110.4 ± 0.3) km/sec and = (2.82 ± 0.1) pc, 
where the uncertainties define the 68% confidence level 
and are derived from the Hessian of the y-square. The re¬ 
sulting black hole mass is M, = (8.0±0.3) 10® Mq and the 
disk mass is approximately equal to the black hole mass. 
We also obtain M = (28.1 ± 0.2) aMQ/yr. If we estimate 
M from the bolometric luminosity (assuming an accretion 
efficiency ry = 0.06), we obtain a ~ 8.3 10“®, a reasonable 
number. 


The best-fit curve resulting from the self-gravitating 


disk m odel is not a power-law. However, Greenhill et al 
( 1996|) were able to obtain a good fit to the data by as¬ 
suming a rotation curve of the form V oc If we 

compute the quantity d InV/d Inr for our best fit model, 
we find that it ranges from —0.35 at the inner edge of the 
disk to —0.30 at the outer edge. 


We have also performed a fit by assuming a Keplerian 
rotation curve. The quality of the fit is definitely worse. 
The minimum y-square in this case is in fact = 1.557 
and is formally rejected at the 95% confidence level. 
The small panel of Fig. ^ compares the two models: the 
Keplerian best-fit curve fails to reproduce both the high¬ 
est and the lowest part of the rotation curve. The resulting 
best-fit value of the black hole mass in the Keplerian fit 
is M, Ri (1.50 ± 0.02) IO^Mq. Note that the total mass 
(disk -I- black hole) of our self-gravitating best fit model is 
roughly the same as the black hole mass of the Keplerian 
fit. Therefore, a non-self-gravitating fit, which attributes 
all the mass to the central object, gives the correct value 
for the total mass, but fails to provide the correct slope of 
the rotation curve. 
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One consequence of the smaller value of the black 
hole mass derived from the self-gravitating disk model 
is that the corresponding Eddington luminosity of the 
black hole is proportionally reduced, hence leading to 
a higher Eddington ratio for the observed central ob¬ 
ject. In fact, based on the results of our ht, we derive 
that iboi/I^Edd = 0.77, to be compared to the value 
T^boi/T^Edd = 0.41 obtained from the Keplerian fit. Note, 
however, that the Eddington luminosity is not well defined 
when the mass is not spherically distributed. 

Some data points (at approximately 0.8 pc from the 
center) are not well fitted by neither the self-gravitating 
nor the Keplerian model. The velocities of these points 
display a higher slope with respect to that predicted by 
the best fit model. This discrepancy could be due to the 
fact that the corresponding masing spots lie on a spiral 
arclet, on a pattern not perfectly perpendicular to our line 
of sight, therefore leading to a smaller projected velocity. 

The statistical significance of the above results, and in 
particular of the uncertainties on the parameters derived 
from the fit, depends on the assumed uncertainty in the 
observed velocities, which may not be easy to estimate. In 
Appendix ^ we discuss in detail this dependence and how 
to discriminate between the different models. 


4.2. Properties of the best-fit disk model 


Eigure g shows the main physical properties of the best- 
fit disk model: the surface density profile, the cumulative 
mass, the equivalent thermal speed, and the aspect ratio 
{h/r), as derived from the self-regulated disk models of 
BL. The cumulative disk mass inside the outer radius of 
the disk is Mdisk(?’out) ~ (8.6 ± 0.6) 10® Mq. Note that at 
large radii tr oc 1/r (also apparent from the approximately 
linear growth of the cumulative mass), with significant de¬ 
viations at r < 0.6 pc. The equivalent thermal speed of 
the self-regulated model is approximately constant at large 
radii, as expected. The number density of H 2 molecules 
in the outer disk is in the range 1 — 5 10®cm“®, com¬ 
patible with the conditions for maser emission and consis¬ 
tent with models of circumnuclear gas heated by the AGN 
( Neufeld et ah , 19^ ; P ier fc Voi^ , |1995| ). X-ray observa¬ 
tions (Bianchi et al, 200l[) an d VLBA radio continuum 
observations (Gallimore et ah, 1997) provide lower limits 
for the electron number density in the disk/torus to be 
rie > 3 10® cm“®, and Ue > 10® ® cm“®, respectively. 


4.3. Consistency of the self-regulated disk assumption 

Starting from Eq. (^ , we can compute the radius at which 
Q is expected to become of the order of unity in a non-self- 
gravitating model. If we refer to the physical parameters 
found from the fit in the previous Section (i.e., M, = 
8 10® Mq, a = 8.3 10“®, and M = 0.23 Mq/jt), we find 
rq 10“® pc. The onset of the self-regulation mechanism 
is therefore expected to take place very deep inside the 
disk, with respect to the radial distances on which we focus 


here, so that the use of a simple completely self-regulated 
model made in the previous Section is justified. 

To strengthen this argument, we have also perfomed 
a fit by assuming the partially self-regulated disk model, 
described in BL, in which the profile of the stability pa¬ 
rameter Q is assumed to decrease according to Eq. (® for 
r rq and to be flat at r rg, with rq = 10“® pc. 
No significant differences are found with respect to the 
completely self-regulated case. 


5. Alternative scenarios 

In this Section we discuss some possible models, alterna¬ 
tive to the self-gravitating disk picture, developed to ex¬ 
plain the non-Keplerian curve in NGG 1068. In particular, 
we will concentrate on the effect of a nuclear stellar cluster 
and on the effect of radiation pressure. Other models also 
exist (for example, the effect of a warp in the disk), but 
have been discussed in the literature less frequently. 


5.1. Nuclear stellar concentration 


One possibility is that the non-point-like source of gravita¬ 
tional field be distributed spherically rather than in a disk. 
A spheroidal nuclear stellar cluster, for example, could 
produce significant changes to the rotation curve traced 
by water maser emission, if its mass enclosed within 1 pc 
from the central engine exceeds IO^Mq ( Kumai| , [1999 ). 
Indeed, a nuclear cusp in the luminosity prohle has 


been observed in NGC 1068 (Thatte et al., 

1997). Based 

on stellar velocity dispersion measurements. 

Thatte et al. 


97D estimate the dynamical mass within 1 (w 69.8 pc 
at a distance of 14.4 Mpc) to be « 6 10®Mq. pchinner^ 
|et al. ( pOOOD report a value of 10® Mq within a 25 pc diam¬ 
eter from the center, based on CO kinematics. It should 
be noted, however, that these mass estimates cannot be 
associated with the stellar cluster only, as the observed 
kinematics includes also the contributions of the central 
black hole and of the disk to the gravitational held. 

Kum^ ( 1999| ) has extrapolated the mass prohle down 
to small radii of the order of 1 pc from the central engine 
of NGC 106 8 starting from the mass estimate of Thatte 
let al. ( 1997 ) and assuming that the stellar cluster can be 
described as a singular isothermal sphere, with stellar den¬ 
sity p oc r“^. In this case, the resulting mass enclosed in 1 
pc would be w 8.6 10® Mq, therefore able to reproduce the 
desired deviations from Keplerian rotation. However, this 
procedure overestimates the stellar contribution in the in¬ 
ner regions, because the stellar cluster is not expected to 
be characterized by the singular isothermal sphere prohle 
down to the innermost regions. Empirically, shallower pro- 


hles are generally found (Faber et ah, 1997). In addition, 


the adiabatic growth of a black hole at the center of a stel¬ 
lar cluster leads to a density prohle p oc ( |Cipollina 

|fc Bertiij , 1994) inside the radius of inhuence of the black 
hole rbh = GM,/a^, where cr* is the stellar velocity dis¬ 
persion (curiously, if we use the value of cr* provided by 
Thatte et al. 1997, we hnd that ri,h is of the order of a 
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Figures. Surface density profile {upper left), cumulative mass {upper right), equivalent thermal speed {lower left), 
and aspect ratio {lower right) of the best-fit accretion disk model to the maser data of NGC 1068. The dotted line in 
the surface density plot is the curve with cr oc 1/r that matches the asymptotic behavior at large radii. 


few parsec, comparable to the scales at which the maser 
emission is observed in NGC 1068). Recent N-body sim¬ 
ulations of the formation of galactic nuclei by merger of 
t wo galaxies with initially steep density profiles {p oc r“^) 
( Milosavljevic fc Merritt , ^001 ) have shown that a shallow 
stellar cusp is left around the nucleus, with p (xr~^ inside 
a break radius of the order of 10^ pc. 

Another difficulty with the picture of a compact stel¬ 
lar cluster is based on dynamical arguments. A collisional 
timescale for a stellar system characterized by a velocity 
dispersion cr* and density p is given by: 


tcoll - 


STrC^pm InA’ 


( 10 ) 


where m is the stellar mass and InA is the Coulomb log¬ 
arithm (we will assume InA = 10). The nuclear stellar 
cluster should have at least a density of p = 10 ^Mq/pc^ 


in its inner regio ns. If we assume cr* = 150 km/sec, the 
value reported by Thatte et al. ( 1997] ), and a stellar mass 
m = Mq, we obtain tcoii ~ 6.5 10^ yrs. Therefore, the 
nuclear stellar cluster, at the high densities required to 
modify the rotation curve, will be subject to collisional 
effects, that may lead to the rapid dynamical evolution of 
the stellar cluster itself. 

To derive firm conclusions about the effect of the stellar 
central concentration on the water maser rotation curve, 
a more detailed knowledge of the properties of the cluster 
at the smallest scales is needed. 


5.2. Radiation pressure support 


Based on radiative transfer models of a thick torus (Pier & 
Krolik, 1992a ^ it has been shown that radiation pressure 
may reduce the importance of the gravitational field of 
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the central object. The relative importance of this effect 
changes at different radii inside the torus and could in 
principle lead to a modification of the rotation curve. The 
main drawback of this picture is that radiation pressure 
is also very efficient at reducing the vertical gravitational 
field. Therefore, we would expect the disk to be rather 
thick, against the observational evidence that the water 
maser emission lies in a geometrically thin structure. In 
fact, the model of Pier & Krolik (1992a) was intended to 
give a theoretical framework for the existence of thick tori 
in AGNs. A more detailed investigation of the effect of 
radiation pressure on thin configurations would therefore 
be needed to assess the importance of this process in the 
case of NGC 1068. 

6. Conclusions 

In this paper we have shown how the study of water maser 
emission in AGNs can be a very useful tool not only to 
estimate the central black hole mass, but also to study 
the properties of the associated accretion disk. In partic¬ 
ular, we have described the non-Keplerian rotation curve 



Figure A. 1. Likelihood ratio A as a function of velocity 
uncertainty AF. The dotted lines indicate the 95% and 
the 99% confidence limits. 


by fitting the data of 

Greenhill & Gwinn 

(1997 

) with the discriminating between the Keplerian disk from the self- 

model described in 

Bertin & Lodato ( 

1999 

). The quality of gravitating disk model. In fact, if the uncertainty in the 


the fit is satisfactory and leads to M, « Mdisk — 8 10® Mq. 
We have also estimated the long-sought value of the a vis¬ 
cosity parameter, obtaining a ~ 10“^. 

Our study has additional interesting consequences. 
Previous estimates of M, , using a non-self-gravitating disk 
model, give M, ~ 1.5 IO^Mq, very close to the value of 
the total black hole -I- disk mass obtained by us. The re¬ 
duced black hole mass obtained in the self-gravitating disk 
scenario leads to a reduced TEdd for NGC 1068, so that 
the Eddington ratio turns out to be larger than previously 
thought. In fact, we obtain Tboi/TEdd — 0.77. More gen¬ 
erally, it is important to derive firm results on M, also 
in view of the correlations recently found between M, 
and the global properties of the host galaxy, such as the 
M, — cr* relation (Ferrarese & Merritt, 2000|; Gebhardt 


As noted at the end of Subsection 4.1, the results of a 


velocities were sufficiently high, we would be unable to 
discriminate between the two different models. 

The self-gravitating disk hypothesis is a generalization 
of the Keplerian one, to which it reduces when M^isk —> 0 
(or equivalently when —>■ oo). Clearly, if the two mod¬ 
els are fitted to the data, the minimum y-square of the 
model with a larger number of parameters is going to be 
smaller. On the other hand, a well-known result of statis¬ 


tical analysis (see |Eadie et al.| |l97l[) enables us to compare 


two competing hypotheses in this case. For Gaussian, in¬ 
dependent measurements, the “likelihood ratio” A, defined 
as: 


A = (minx^)Kep - (minx^)s 


(A.l) 


et ah, 2000a). In this context, it is worth noting that the 
masses derived from water maser emission have been con¬ 
sidered as the most reliab le estimates; NGC 1068 is indeed 
one of the galaxies in the Gebhardt et al. (2000a) sample. 
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Appendix A: A likelihood ratio test for the 
self-gravitating disk model 


statistical analysis of the water maser data depend signif¬ 
icantly on the uncertainty assigned to the velocity data 
points. In this Appendix we examine this issue and dis¬ 
cuss how it is connected with the general problem of 


is distributed like a x-square with n degrees of freedom, 
where n is the number of new parameters in the more gen¬ 
eral model (in our case n = I), under the hypothesis that 
the less general model is correct (in our case, the Keplerian 
model). Of course, the resulting value of A depends on the 
assumed uncertainties. Therefore, an interesting question 
is how large should the velocity uncertainty (for example, 
that associated with the magnitude of turbulent veloci¬ 
ties) be to make the Keplerian hypothesis acceptable with 
respect to the self-gravitating one. 

In Fig. A.l the likelihood ratio A is plotted as a func¬ 
tion of the assumed velocity uncertainty AV. The two dot¬ 
ted lines define the 99% and 95% confidence limits. The 
Figure clearly shows that the Keplerian model should be 
rejected against the self-gravitating one, with 99% con¬ 
fidence, if AV < 25 km/sec, and with 95% confidence, 
if AV < 30 km/sec. Note that the required uncertainty 
to make the Keplerian model acceptable is more than 20 
times the formal instrumental uncertainty. This may be 
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the reason why Greenhill & Gwinn (1997) argued that 
“the scatter in the data may indicate turbulent velocities 
of up to a few tens of km/sec”. 


References 

Baan, W. A. & Hashhick, A. 1996, ApJ, 269, 473 
Bardou, A., Heyvaerts, J., & Duschl, W. 1998, A&A, 337, 
966 

Bertin, G. & Lin, C. 1996, Spiral Structure in Galaxies: a 
Density Wave Theory (Cambridge: MIT Press) 

Bertin, G. & Lodato, G. 1999, A&A, 350, 694, (BL) 
Bianchi, S., Matt, G., & Iwasawa, K. 2001, MNRAS, 322, 
669 

Cipollina, M. & Bertin, G. 1994, A&A, 288, 43 
Eadie, W. T. et al. 1971, Statistical Methods in 
Experimental Physics (London: North Holland) 

Faber, S. M. et al. 1997, AJ, 114, 1771 
Ferrarese, L. & Merritt, D. 2000, ApJ, 539, 9L 
Gallimore, J. et al. 1997, Nature, 388, 852 
—. 2001, ApJ, 556, 694 
Gebhardt, K. et al. 2000a, ApJ, 539, 13L 
—. 2000b,AJ, 119, 1157 

Greenhill, L. & Gwinn, C. 1997, Ap&SS, 248, 261 
Greenhill, L. et al. 1996, ApJ, 472, 21L 
Hure, J. 2000, A&A, 358, 378 
Kumar, P. 1999, ApJ, 519, 599 
Lin, D. & Pringle, J. 1987, MNRAS, 225, 607 
Magorrian, J. et al. 1998, AJ, 115, 2285 
Milosavljevic, M. & Merritt, D. 2001, ApJ, 563, 34 
Miyoshi, M. et al. 1995, Nature, 373, 127 
Neufeld, D. A., Maloney, P. R., & Conger, S. 1994, ApJ, 
436, L127 

Paczyhski, B. 1978, Acta Astronomica, 28, 91 
Pier, E. & Krolik, J. 1992a, ApJ, 399, 23L 
—. 1992b, ApJ, 401, 99 
Pier, E. & Voit, G. M. 1995, ApJ, 450, 628 
Pier, E. et al. 1994, ApJ, 428, 124 
Pringle, J. 1997, MNRAS, 292, 136 
Schinnerer, E. et al. 2000, ApJ, 533, 850 
Schmitt, H. et al. 2002, ApJ, 575, 150 
Shakura, N. & Sunyaev, R. 1973, A&A, 24, 337 
Thatte, N. et al. 1997, ApJ, 490, 238 
van der Marel, R. P. & van den Bosch, F. 1998, AJ, 116, 
2220 





